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Abstract 



We derive an analog of the master equation, obtained recently for correlation 
functions of the XXZ chain, for a wide class of quantum integrable systems 
described by the i?-matrix of the six- vertex model, including in particular con- 
tinuum models. This generalized master equation allows us to obtain multiple 
integral representations for the correlation functions of these models. We apply 
this method to derive the density-density correlation functions of the quantum 
non-linear Schrodinger model. 



^LPTM, UMR 8089 du CNRS, Universite de Cergy-Pontoise, France, kitanine@ptm.u-cergy.fr 

Laboratoire de Physique, ENS Lyon, France, karol.kozlowski@ens-Iyon.fr 
^ Laboratoire de Physique, UMR 5672 du CNRS, ENS Lyon, France, maillet@ens-lyon.fr 
* Steklov Mathematical Institute, Moscow, Russia, nslavnov@mi.ras.ru 

^ Laboratoire de Physique, UMR 5672 du CNRS, ENS Lyon, France, veronique.terras@ens-lyon.fr, on leave 
of absence from LPTA, UMR 5207 du CNRS, Universite Montpellier II 



1 



1 Introduction 



The method to study the spectrum of the Heisenberg antiferromagnet was proposed by H. 
Bethe in 1931 [Tj. This method, called later Bethe ansatz, was found to be very effective for 
a complete description of the spectrum of quantum solvable models. During the ensuing years 
the Bethe ansatz was successfully extended and applied to the study of various problems (see 

e.g. [21 El H El El E]). 

A further development of this method was done in [S], where the algebraic version of the 
Bethe ansatz was formulated. In the algebraic Bethe ansatz one deals with a set of operators 
(the operators entries of the quantum monodromy matrix) possessing quadratic commutation 
relations (governed by the i2-matrix solving the Yang-Baxter equation). Different representa- 
tions of this algebra correspond to different physical models associated to the same i2-matrix. 
Thus, the algebraic Bethe ansatz allows one to consider a wide class of quantum integrable 
models from a uniform algebraic viewpoint (see [9] and references therein). 

In spite of its successful application to the calculation of the spectrum of quantum Hamilto- 
nians, the possibility to use the Bethe ansatz to evaluate the corresponding correlation functions 
remained involved for a long time. Basically the first complete results for integrable models cor- 
relation functions concerned the models equivalent to free fermions [SI [TU\ \TI\ [T^ [THl [2] for 
which different methods apply. Attempts to compute the correlation functions of local opera- 
tors by the algebraic Bethe ansatz for more general models were undertaken in [151 IlSl El E] • 
One of the main problems there is to embed in an appropriate way the local operators into the 
quadratic Yang-Baxter algebra describing the spectrum of these models. 

The solution of the quantum inverse scattering problem, first obtained for the XX Z chain 
|19j . leads to a neat solution to this question: the local operators are then given as explicit 
products of the quantum monodromy matrix entries. Hence it gives a powerful tool for the 
computation of the correlation functions via the algebraic Bethe ansatz (see e.g. |20j and further 
developments reviewed in [H]). In turn, the explicit and compact expressions for local operators 
in terms of the entries of the quantum monodromy matrix reduce finally the calculation of the 
correlation functions to the well studied problem of the computation of scalar products in the 
algebraic Bethe ansatz framework [26 1 127 1 [T9j. It is worth mentioning, however, that the solution 
of the quantum inverse scattering problem has been obtained so far only for lattice models [22], 
the case of continuum (field theory) models remaining an open question. 

The main goal of this paper is to extend some results obtained for the XX Z chain to a wider 
class of integrable models described by the algebraic Bethe ansatz, but for which the explicit 
solution of the quantum inverse scattering problem is not known at present. Namely, we derive 
an analog of the master equation for correlation functions of the XX Z chain obtained in [23] 
for the generalized model proposed in [241 [T5l I16j . This model includes as a particular case 
continuum (field theory) models, the quantum Non Linear Schrodinger (NLS) model being one 
typical example. Just like in the case of XX Z chain, this master equation can be used to derive 
multiple integral representations for a certain class of correlation functions. In particular, we 
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obtain the density-density correlation functions of the quantum NLS model. The basic idea 
to reach this goal is to find a way to embed certain local operators of these models into the 
Yang-Baxter algebra of the quantum monodromy matrix entries, or at least, into the algebra 
of the quantum monodromy matrix entries corresponding to a subpart of the model. This can 
be achieved in particular for the number of particle operator associated to a segment [0, x\ 
in the quantum NLS model. Using then algebraic Bethe ansatz techniques, together with the 
determinant formula for the scalar products of Bethe states |26[ l27j. it is possible to compute the 
corresponding correlation functions in a rather compact way, hence leading to density-density 
correlation functions. The same technique can be applied for more general correlation functions 
as well. 

The content of the paper is as follows. In Section 2 we recall some basic definitions and 
formulas necessary to study the generalized model in the framework of algebraic Bethe ansatz. 
In Section 3 we consider the two-site generalized model and derive the master equation for 
the expectation value of the operator of number of particles in the first site. In Section 4 we 
use this master equation to evaluate the multiple integral representation for the density-density 
correlation function of the Quantum Nonlinear Schrodinger model. The proofs of some algebraic 
identities are collected in three appendices. 



2 Generalized model 



In this section we recall briefly the main definitions of the algebraic Bethe ansatz and introduce 
necessary notations. We refer the reader to [8], [9] for more details. 

The main objects of the algebraic Bethe ansatz are the i?-matrix R{X), the monodromy ma- 
trix T(A) and the pseudovacuum vector ]0). We consider quantum integrable systems associated 
to the ii!-matrix of the six- vertex modejj 



RiX) 



( sinh(A -|- vi) 



V 





sinh A 
sinhr/ 






sinh 7/ 
sinh A 











(2.11 



sinh(A + rj) J 



where A and r/ are complex numbers. Usually the parameter r] is related to the coupling constant 
of the Hamiltonian. 

The monodromy matrix of such model is a 2 x 2 matrix 



T(A) 



^(A) B{X) 
C(A) D{X) 



(2.2) 



with operator-valued entries A,B,C and D which depend on a complex parameter A and act 
in the quantum space of states of the Hamiltonian. These operators satisfy a set of quadratic 



^One can also consider the particular case where hyperbolic functions in the i?- matrix degenerate into the 
rational ones. 
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commutation relations given by 

i?(A - M)nA)r(^) = Tip)T{X)R{X - /.). (2.3) 

Here r(A) = T(A) / and T(A) = / O T(A). 

The pseudovacuum |0) is characterized by the action of the entries of the monodromy matrix 
on this vector 

A(A)|0) = a(A)|0), D(A)|0)=d(A)|0), C(A)|0)=0. (2.4) 

Here a(A) and d{X) are complex functions depending on the specific model. In the framework 
of the generalized model they remain free functional parameters. Various models associated to 
the same i?-matrix are described by different representations of the ratio r(A) = a{X)/d{X). 
The dual pseudovacuum (0| = |0)^, (0|0) = 1 is characterized by similar equations 

(0|.4(A) = a(A)(0|, {0\D{X) = d{X){0\, {0\B{X)=0. (2.5) 

The trace of the monodromy matrix T^/j,) = A(/x) + D{fj,) appears to be the generating 
function of the integrals of motion of the model. In particular the eigenstates of T(iJ,) coincide 
with the eigenstates of the Hamiltonian. We shall also consider a more general object 

T,(/i) = A{f,) + ^D{f,), (2.6) 

which is called the twisted transfer matrix. Here k is a complex parameter. The eigenstates 
(and their dual ones) of the operator 7^(/i) can be written in the form 

N N 

IMW)) = n ^(^^OIO), (V'.({A})| = (0| n C{X,), N = 0,1,... (2.7) 

i=i i=i 

where the parameters Ai, . . . , Xn satisfy the system of twisted Bethe equations 

J,(A,|{A}) = 0, j = l,...,N. (2.8) 

Here 

N N 

3^«(/«|{A}) = a{fi) n sinli(Afc - /i + r?) + Kd{fi) sinh(Afc - /i - r?). (2.9) 

k=l k=l 

The corresponding eigenvalue of Xt(/") on \tpK{{^})) (or on a dual eigenstate) reads 

N . , , . N 



^J-^ smh(Afc-//) smh(//-Afc) 

Formulas (|2.7p - (|2.10p remain valid for the particular case k = 1 for which we omit the subscript 
K, denoting for example, \^p{{X})), 3^(/u|{A}) and t(;u|{A}). 
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A central role in the calculation of the correlation functions is played by scalar products of 
states 

TV N 

SNmm) = {o\i{c{f^,)i[B{x,m, (2.11) 

where the parameters /ii, . . . , /iat and Ai, . . . , Aat are arbitrary complex numbers and generically 
do not satisfy any constrains. Using the commutation relations (|2.3|) one can obtain an explicit, 
but not very compact formula for this quantity |15|, [U] ■ 

(01 fr fr :^|o) = fr I 

y difij) y diXj) ' ' ii sinh(A,- - Afc) sinh(/xfc - 

X y^[-i)iPi-)\+[ni)\ rr ^(^j rr det t(A,-, /x,) det a,) 

X n n ^(^^" n n ^(-"^^ n n ^^^^^ k) n n ^(^^^^ (2-12) 

Let us describe the notations used in (j2.12p . First, we have introduced the two functions 

Slllh TJ 

h{X,ii) =smh{X- 11 + rj), t(A, //) = —— . . , ■ — r- (2.13) 

smh(A — /i) smh(A — fi + 7]) 

The set of the parameters {A} is divided into two disjoint subsets {A} = {A}q U {X}a- Similarly 
{fi} = {/ij-y U {/J.}^. Hereby #{X}a = #{^}^. The parameters in each subset are ordered in 
the natural way, namely {Aq^ , Aqj, . . . } if ai < 02 < . . . etc. The sum in (|2.12p is taken with 
respect to all possible partitions of this kind. The symbol P{a) (respectively P{'y)) means the 
permutation P{{a}, {a}) = 1, . . . ,N (respectively ^({7}, {7}) = 1, . . . , A^). Similar notations 
will be used through all the paper. 

The formula (I2.12|) can be drastically simplified if one of the states is an eigenstate of the 
twisted transfer-matrix. Let us introduce, for arbitrary positive integers n and n' (n < re') 
and arbitrary sets of variables Ai, . . . , A„, fii, . . . , fin and z^i, . . . , Vn'i such that {A} C {z/}, the 
following n X n matrix ^^({A}, {/ulll^^}) 



a=l 



(f^K)jfc({A}, {/i}|{z^}) = a{fjLk)t{Xj,fik) Y\^smh{va - fik + rj) 

n' 

:d{fik)t{fik,Xj) JJ sinh(i^a - /"fc - ??)• (2.14) 



a=l 
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Then, if {A} satisfies the system of twisted Bethe equations, 



N N 

{0\llC{^^,)\M{m = (VK({A})in^(/^.)|0) (2.15) 
i=i j=i 

- ,N . . ' . \ . , ^-detf^.aAI.MKA}). (2.16) 



nj>fc sinh(Aj - Afe) sinh(/Xfc - [ij] 



N 



Different proofs of this formula for the case k = 1 were given in [26^ [27] , and for the particular 
case of the XX Z chain in [19j. For completeness we give in Appendix A another proof for 
arbitrary k, which will be used later. 

Observe that formula (j2.15p implies two different representations for the particular case 
of the scalar product {tP{{ii})\iPk{{^})) , for which in addition {fi} is a solution of ordinary 
(non-twisted) Bethe equations: 

= . . ^ ■detl7«({A},M|{A}) (2.17) 



n,>fc sinh(Aj - Afc) sinh(/ifc - fij) 



N 



'■J 



T-lN 1/ 

N . . \, :-detJ7(M,{A}|M), (2.18) 



n,>fc sinh(Aj - Afc) sinh(/XA; - /Uj) 



N 



where VL = ^J^^^^. Indeed, one can consider either the state {i^{{^J'})\ or the state |?/'k({A})) as 
a particular case of an arbitrary state. We conclude that 

N N 

n d{Xj) ■ det n^{{X}, {^}|{A}) = n din) ' det 0({^}, {A}[{/i}). (2.19) 

This identity holds if the parameters {A} (respectively {/x}) are solutions of twisted Bethe 
equations <\2.8h (respectively ordinary Bethe equations with k = 1), see Appendix[B]for a direct 
proof. Otherwise both sides of the equation (j2.19p define different functions of {A} and {/x}. 
Setting fij = Xj, j = 1, . . . , N in (j2.15p one has 

{MimMim = ^^^=''^^^^^ ■ det (^-A-y,{X^\{X})^ , (2.20) 

n sinh(Aa - Afe) 

a, 6=1 

where is given by ()2.9p . 



3 Two-site model 

A number of correlation functions can be computed in the framework of the generalized two-site 
model introduced in |15i I16j (see also ^). In the two-site model the monodromy matrix T(A) 
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is presented as the product of two operators T(A) = T2(A)ri(A), where 



The commutation relations between the entries of each Tj{\) are given by the i?-matrix (j2.ip . 
the operators from the different sites commuting with each other. The monodromy matrices 
Tj(X) have their own pseudovacuum vector |0)j and dual one Hereby |0) = |0)2 (8) |0)i. The 
action of the operators on these vectors is similar to (|2.4p . (j2.5p 

Aj{X)\0)j = a,-(A)|0)„ Dj{X)\0) = d,(A)|0),-, C,(A)|0),- = 0, 

(0|,.4,(A) = a,(A)(0|„ {0\jD,{\) = d,(A)(0|„ (0|,S,(A) = 0. 

It is easy to see that a(A) = ai(A)a2(A) and d{X) = di{X)d2{X). Below we shall also use the 
functions 

Following the work [15] we introduce now the operator Qi of number of particles in the first 
site. This operator is defined by its commutation relations with the entries of the monodromy 
matrices Tj{X) 

[Qi, Bi{X)] = Bi{X), [Qi, Ci{X)] = -Ci{X). (3.3) 
All other commutation relations are trivial. The action of Qi on the pseudovacuum is given by 

Qi|0)i = 0, (0|iQi = 0. (3.4) 

The expectation value of the operator e^^^ , where /? is a complex number, is a generating 
function for correlation functions. In order to compute this expectation value in the framework 
of the generalized two-site model one should use the commutation relations (13. 3p . We use the 
result obtained in [151 [E] • Let 

m)) = (oi n TTTT • • n tMio) • n - ^>^) ^^^^('^'^ - /^i)' ^ = (^.s) 
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where {A} and {//} are generic complex numbers. Then 



X sinh(Aa - Af,) 

a>b 



sinh(/Ua - Hb) 

a<h 



sinh(Aa - Xb) 

a>b 



sinh(/ia - fib) 

a<b 



ago ae7 aG« fego 067 6G7 

Here n denotes the number of elements in the subset {a} (or {7}). Thus, we see that this 
matrix element can be reduced to the scalar products in the first and in the second sites. For 
both of them one can use the representation (|2.12|) with appropriate replacement of r(A) by 
/(A) or m(A). 

Substituting (I2.12|) for each scalar product in ()3.6p we obtain a sum over 4 partitions of the 
set {A} and 4 partitions of the set {^}. It was shown in the papers [151 [E] that the sum over 
four of these partitions can be calculated, if the sets of the parameters {A} and {fj,} are two 
different solutions of Bethe equations. This result, however, does not allow one to set {A} = {fj,} 
in this simplified formula, since two different solutions of Bethe equations can not be sent to 
each other. We consider therefore a more general case, when only the set of the parameters {A} 
satisfies Bethe equations, while the parameters {/x} remain arbitrary complex numbers. In a 
sense this case is analogous to the one considered in Section [2] for scalar products. 

The main result of this section and of the paper is the generalized master equation for the 
matrix element of the operator e^*^^. 

Theorem 3.1. Let the parameters Ai, . . . , Xn satisfy Bethe equations at k = 1, while fii, ■ ■ ■ , fiN 
are generic complex numbers. Then 



r(Mu{A}) 

detjv n,{{z}, {fi}\{z}) detN n^{{z}, {X}\{z}) 

where the integration contour T ({fj,} L) {X}) surrounds the poles of 0,^ in the points {A} and {/x} 
and no other pole of the integrand. 
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The proof is given in Appendix O We would like to stress here that the only information 
used to derive this formula are the commutation relations (I3.3p . 

Note. In the XXZ model the two "sites" of the generalized model correspond to the 
division of the chain into two subchains. In this framework, the operator e^'^^ has an explicit 
representation in terms of local spin operators, 

^,Q. = n + ^ • <) = ff - ' ^'=-^ (3.8) 

n=l ^ ^ \ n=l / 

where m is the length of the first subchain. Using the solution of the inverse scattering problem 
pO] one can obtain the equation (13. 7p by the method described in p^. However, this approach 
fails in the case of the generalized model, since the solution of the inverse scattering problem 
for such a model is not known at present. 

Since the parameters {fi} in the equation (j3.7p do not satisfy any constraints one can consider 
the particular case = Xj, j = 1, . . . ,N. In this case we obtain for the normalized expectation 
value of the operator e^^^ 

r({A}) 



det7vf^«({4,{A}|{z}) 



n3^.(z,|{z})det^^^%P 



(3.9) 



4 Quantum Nonlinear Schrodinger equation 

As an application of the master equation ()3.9p we derive now a multiple integral representation 
for the density-density correlation function of the one-dimensional Bose gas. The method de- 
scribed below is very similar to the one used in [231 [29] for the XXZ chain, therefore we refer 
the reader to the mentioned papers for more details. 

The model of the one-dimensional Bose-gas is described by the Quantum Nonlinear Schrodinger 
equation. The Hamiltonian of this model on the ring of length L is given by 

H = J (j3x'ipH^)9xi^{x) + c^l^'^ {x)^p\x)^p{x)^p{x) — h^\x)'ilj{x)j dx, (4.1) 

where c > is the coupling constant and h > the chemical potential. The operators tp'^ and i/j 
are canonical Bose-fields obeying standard equal-time commutation relations. This model was 
studied via the algebraic Bethe ansatz in PH ISl M- The operator Qi giving the number of 
particle in the interval [0, x] for this model has the following explicit representation 

Qi= [ i'\z)iP{z)dz. (4.2) 
Jo 
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The commutation relations (13. 3p can be proven directly. Then, the ground state expectation 
value of e^^^ can be used to compute the density-density correlation function [15| 

(^t(^)^(^) .^t (0)^,(0)) = \^ ^(e/^Qi)|^^^. (4.3) 

The master equation (j3.9p permits us to obtain in a simple way a multiple integral representation 
for this quantity in the thermodynamic limit. 

In the case of the Quantum Nonlinear Schrodinger equation one has 

a(A) = e-'^^l^, d{\) = e'^^l'^, = e-^^\ (4.4) 

One should also replace all hyperbolic functions in 0^ and by the rational ones and set 
r] = —ic. The system of Bethe equations then reads 

\j -Xa-lC 

a — 1 

and all its solutions are real numbers. For the solution describing the ground state of the model 
the spectral parameters {A} belong to an interval [q-,—q]i where q depends on the coupling 
constant c and the chemical potential h. The system of twisted Bethe equations can be written 
in the form 



n ^' , ^ = e^- (4.6) 



a = l 



Therefore all the roots of the twisted Bethe equations belong to the line M + i(3/L. 
The equation ()3.9p takes the form 



«n - (zlV: i n (i?a .,<,.,.-.„.,-.,,\ , ^[<le%».({4.{A>l{.-})]' ^^^^^ 



r({A}) Y{yK{zj\{z})detN 



dy{\,\{\}) 



dXi, 

Let /3 < 0. We will now transform this representation by several moves of the integration 
contour. Note that the integrand vanishes when \zj\ — > oo. Therefore, one can evaluate the 
integral ()4.7p by the residues at the poles of the integrand outside the integration contour. 
These poles correspond to the solutions of the twisted Bethe equations (14. 6p . which lie on the 
line M + i/3/L. Thus, as soon as we have moved the contour to these roots, we can replace one 
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of det by det O due to (I2.19p , and we have 
V, i - - \ 



<^^*>=Hn / - 

\jR+i/3/L-iO R+j/3/L+iO / 



N 



-ix{zj—Xj) 



det^ a^{{z}, {\}\{z}) det^ 17({A}, {z}\{\}) ^ ^^ ^^ 
n3^«(z,|{z})det^^^%P 

The integrals over M + i/3/L — iO vanish due to the exponential decreasing of the factors 
exp(— izjx) in the lower half plane. Thus, we obtain 



(-1) 



N 



N , N 



-ix{zj 



N\ J 27ri 

R+iP/L+iO 



detM n^{{z}, {\}\{z}) detN {z}\{\}) ^ 



N 



n3^«(^,[M)det^^ 



9y(A,l{A}) 

k 



It remains to move the contour up, through the real axis. Note that the matrix r2({A}, {-z}|{A}) 
has no poles at Zj = Xk, therefore we cross only the simple poles of detAr risdz}, {A}|{z}). 
Taking into account that 



Res det nMz},{\}\{z}) 

N 



{-ir+'y.{Xb\{zy^z,},x,) 



X det / z4, {A ^ Xb}\{z ^ z4, Afe), (4.10) 



we arrive at the sum over partitions 



^ (_1)^ 



n=0 



dz 



27rz 



^ det„ Vtn{{z}, {X}^\{z} U {A}^,) detjv Vt{{X}, {z} U {A}^,|{A}) 



n3'.(^,|{4u{A}e.)det^™ 



gy(A,i{A}) 



The representation (j4.1ip is the complete analog of the one obtained for the XXZ chain (see, 
for example, [201 [281 [231 [29] ) . Therefore the remaining steps of the computation, in particular 
the thermodynamic limit, are the same as the one's in the XXZ model. We refer the reader to 



11 



these cited papers for the details and give here the final result for the ground state expectation 
value of the operator e^^^ in the thermodynamic hmit 



q 

OO ^ J. r ''^ rl " 



1 f f n ^^J TT (Ag - Zb- ic){zb - Xg- ic) 

^. 1 (-^a - - ic){Za - Zb- ic) 
a,b=l 



° -q R+iO 

n 

^TJ^-^xiz,-x,) detM({A},{4).det[p(A„z,)]- (4.12) 

6=1 

Here 

%({A}.M)-*fa.A,)+e^*(A,..on ^;;;;;^:;^;^;;°; . (4.13) 

and the function p(A, z) solves the integral equation 

- y (A - ^)^ + = (A-z)(A-.-zc) - ('-^'^ 

Observe that the series (|4.12p is completely analogous to the one obtained for the XX Z 
chain in [28j. The function p{X, z) is just the rational limit of the corresponding quantity for 
the XX Z model. It is a universal function depending essentially on the i?-matrix. 

The series (|4.12|) is especially effective in the large coupling limit c — > oo. Indeed, one can 
see that only the first m + 2 terms of the above series contribute to the m — th order in 1/c- 
expansion of the density-density correlation function. In particular it reproduces the known 
results at the free fermion point c = oo (see e.g. ^ and references therein). 



Conclusion 

We have shown that the analog of the master equation for correlation functions exists for a wide 
class of integrable systems described by the i?-matrix of the six-vertex model. Our method 
applies, in particular, to the Quantum Nonlinear Schrodinger equation and the lattice Sine- 
Gordon model. Like in the case of the XX Z chain, this generalized master equation provides 
the direct analytical link between multiple integral representations for the correlation functions 
and their expansions over form factors. Using the generalized master equation one can also 
obtain various integral formulas for the correlation functions in the thermodynamic limit. In 
the forthcoming publications we are planning to give such integral representations, which are 
convenient for the asymptotic analysis of correlation functions. 
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A Proof of the determinant representation for the scalar pro- 
duct 

Consider an auxiliary function G^^^ depending for given n and N {0 < n < N) on three sets of 
complex parameters {^i, . . . , {i^i, . . . , vw-n} £^nd {Ai, . . . , A^v} 

GPm,M,W)= E (-1)'''^"^' det i(A,-,a)- det iK,A,) 

«U<5 fc=l,...,n fc=l,...,Ar-Ti 

#{a}=n 



X 



' n N—n n N—n 

n n ^(^^" n n ^(^b, - n n ^(^^ > n n ^(^»' 

^a=l 6=1 aeabea a=lbea a=l bea 



n N—n 

X n n ^(^'" n n ^(^«' )■ (^-i) 

o=l bea o=l bea 

Here the sum is taken with respect to all partitions of the set {A} into two disjoint subsets 
{A} = {A}q; U {X}a, such that #{«} = n. 

Lemma A.l. For arbitrary complex numbers {^i, . . . , {vi, . . . , fN-n} 0,'f'd {Ai, . . . , Aat} 

G&^({e},M,{A})^0. (A.2) 

Proof. For N = 1 the statement of the lemma is evident, since either n = OorA?^ — n = 0. 
Let G^^^^d^}, {u}, {A}) = 0. Consider the properties of G^^^ as a function of Ai, . . . , Ajv, with 
other parameters fixed. Obviously G^^^ = as soon as Aj = Afe {j,k = 1, . . . , N, j ^ k). Hence, 
the combination 

G'^>m,M,m = n sinh(A, - a,) ^^-^^ 

k<j 

~ (n) 

has no poles at Aj = A^. Therefore the only possible singularities of G)^^ are the poles at Xj = 

~ (n) 

and Xj = Uk- It is easy to see that the residues in these poles are proportional to G]y^_^, for 
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instance, 



ResGi;)({e},M,{A}) =G);'i,({e/ei},M,{A/Ai}) 

Al=Cl 



~ (n) 

Due to the induction assumption these residues vanish. Thus, G)^ is a holomorphic function 
of each Xj in whole complex plane. It remains to observe that {{(,}, {u} , {X}) IS an ZTT- 
anti-periodic function of the parameters {A} and that it exponentially decreases at |KAj| —>■ oo. 
Hence, GS;?)({e},M,{A}) =0. □ 

Consider now the equation (|2.12p . Let one of the states, for instance Ylf=i B{Xj)\0), be an 
eigenstate of the twisted transfer-matrix. Then, due to the twisted Bethe equations (12. 8p . one 
can express r{Xj) as 

Substituting this into (j2.12p we obtain 

N , N s N 



\ ^if^j) fj^ d{Xj) j-J-^ sinh(Aj - Afc) sinh(/ifc - ijj) 

X V(-1)[^(")] + [^W] (^-l)nN-n^N-n TT ^(^^^ ^ ) 

7U7 

^ n n n n n n ^^^-^ n n ^^^b, /^a), (a.6) 

where n = #{a}. For given partitions {7} and {7} one can set {/ij-y = {^^i, . . . ,Cn}) {lJ'}'y = 
{ui, . . . , z/jv_„} and then apply Lemma [A. II 



(»in§^n§Mi«>=n 



_i JJ- d(Aj) \\ sinh(Aj - Afc) sinh(^fc - /ij) 

X V(-1)I^(-)]+I^W](-1)"^-"k^-" rr r(^,) dett(A„/ifc) dett(/xfc. A,-) 

7U7 

N N 

xl{llh{Xb,fia)l{llh{^,b,Xa)- (A.7) 

0676=1 a=l 6e7 
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It remains to use the Laplace formula for the determinant of the sum of two matrices and we 
arrive at ()2.15p . 



B Identity for determinants 

Let 

Gjk = t{nk, Xj) Y[ h{fJ-k,Xa) + Kt{Xj,flk) Y[ h{Xa,fJ.k) JJ ^' (B.l) 
a=l a=l a=l V^^'^^'' 

N N ^ h(\ X ) 

Gjk = tifik, Xj) n Hfia,Xj) + Kt(A„ ^fc) n h{X,,f,a) n , /\ (B.2) 

a=l a=l a=l ■' 

Lemma B.l. For arbitrary complex numbers Ai, . . . , Xn and ni, . . . the following identity 
holds 

det G = det G. (B.3) 

N N 

Proof. The proof of this lemma is similar to the one of Lemma [A. 11 For = 1 the statement 
of the lemma is correct. Assume that it is still correct for A^ — 1. Then, considering the equation 
(|B.3p in the points Xj = /ifc, we see that the corresponding residues are proportional to the 
determinants of the size A^ — 1, which are equal to each other due to the inductive assumption. 

The only nontrivial step is to prove that the determinants of G and G have no other singu- 
larities except at the points Xj = ^k- At the first sight one could expect that detG might have 
simple poles at h{Xj,Xk) = 0, and det G at h[^j,fik) = 0. It is easy to see, however, that the 
corresponding residues vanish. Indeed, let, for example, /i(Ai,A2) = 0. Then the first line of 
the matrix G becomes singular 

^ ^ h(\ \ \ 

Gik ^ ^t{Xi,f^k) n ^(Ai,/xa) n TTtVt- 

a=l a=l 

On the other hand one has for the second line of G 

N 

Gsfc ~t(/ifc,A2)n^(Ma'^2). (B.5) 

a=l 

Taking into account that t(/iA;,A2) = t(Ai,;Ufc) at h{Xi,X2) = 0, we see that the first and the 
second lines of G are proportional to each other. Hence, the corresponding residue vanishes. 
Similarly one can prove that detG has no poles at h[^j^^k) =0. □ 

The identity (jRT]) allows one to prove directly the equation (j2.19p if the parameters {A} 
(respectively {^}) are solutions of twisted Bethe equations (12. 8p (respectively ordinary Bethe 
equations with k = 1). 
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C Proof of the theorem 13.1 

We start with the equation (13. 6p . Substituting (12.120 for each scalar product we obtain a sum 
over 4 partitions of the set {A} and 4 partitions of the set {/i}: 

^(^?)) = E n ^(^") n "^(^") n ^(^«) n ^(^'^) 

7^ U7_ U7_ 

^ ^n++?i_ t{Xj,fik) det t{fik,)^i) det t{Xj,^k) det t[fik,Xj) 

fce=7-^ fceE7_ fc£7-{- fce7_ 

X J] J] /i(Ab,/za) n n ^('";"^«) n n ^^^^ao) n n ^(z";"/"-) 
^ n n ^(^'"/"«) n n ^^i^^Aa) n n MAfe,Aa) n n ^i^b^i^a) 

aecf+Ua. feea+Ua_ aG7+U7_ be7+U7- 

Here, as usual, P{a) is the permutation acting as P({a+}, {«-}, {«+}, {«-}) = 1, • • • , and 
the permutation P{'^) acts in a similar way. The numbers of elements in the subsets {«+} and 
{q;_} are denoted as n+ and n_ respectively. 

The sets {^u} and {A} in (fUl^ are still generic complex numbers. If the set {A} satisfies 
Bethe equations (|2.8|) with k = 1, then we can express r[\a) similarly to (jA.Sp . Denoting 

ao = a+ U q;_, = a+ U a_, 

70 = 7+ U 7_, 70 = 7+ U 7_, 
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we obtain 



U7_ U7^U7_ 



X det t{Xj,^k) det t{fik,Xj) JJ JJ h{Xa,\b) JJ JJ Hl^b,l^a) JJ h{na,fJ'b) 

ke-y^ fcS7_ aeSo bGao be7- aG70 ''€7+ "G70 



X n J] h{Xh,fia) n n ^(^^f'^'^) n n ^iXaAb) n n ^(^-z^^') 

X TT det t(/^fc,A,) det t{X,,^ik) 

X n n ^(^f'^'^) n n ^^^f^a) n n ^(^a'^^-) n n Hpa,f^b)- (c.3) 

aeo- be-y- 1G7+ feG(5+ ago- feG5:+ aG-y- feG7_|. 

The sets {A}q,_ and {A}a_^ enter separately only the two last lines of (jC.Sp . Hence, we can sum 
up over these partitions (at others partitions fixed) and apply Lemma lA.ll 



y det t{fXk,Xj) det t{X,,fik) 

0_Ua+ fcS7_ fc67+ 



X n n ^(/^'"^«) n n h{h,^^a) n n ^(^^'^t) n n h{i^a,f^b) 

aea- b€y- a-&l+ bea+ aG«- bGa+ aG7_ 6e7+ 



V det t(/Xfe,A,) det t(A„/Xfc) 

Q_Uq+ fee7_ fc67+ 



X n n K^-^^'i>) n n ^(/^^"^-)- (c-4) 

fe€7+ a€i5o 6G7_ aSao 



Substituting this into (|C.3p we see that the sum over the partitions {q_}, {«+}, {7-}, {7+} 



17 



jives exactly the Laplace decomposition of the determinant. Thus, we obtain 



m)) 



nKfJ-a) IT "i(-^a) • /^^^ • det t{Xj,iJ,k) det t{nk,\j 



7+U7-U70 



fee 7-^ 



det 



fee7o 



aG7o 



ae7o 



n n ^(^-^f) 



Now, for arbitrary fixed partitions qq and 70, we can smn up with respect to the remaining 
partitions. Using again Lemma lA.ll we have 



V . det t(A„/ifc) det t{i,k,Xj)- n n 



a_j_ UcK_ 
7+U7- 



feg7_|_ 



fce7- 



ae7+ {,£70 



X n n hih,f^a) n n ^(z^^"^-^) n n hi\a,Xb) n n ^(/^^"/^») 



V . det t{Xj,i,k) det t(/i,,A,)- TT TT 



U<:k_ 
7+U7_ 



fe€7+ 



h{fla,flb) 



X n n ^(^"'/^fc) n n Ki^hAa) 

be7+ aGao 6G7_ aGoo 



det 

*:e70 



K^^a,|J'k) 



(C.6) 



18 



Applying now Lemma IB. II we finally arrive at 

m)) = E n n -(^'^) • n n ma., a,) 

a€7o aGoo aGao feGoo 



70^70 



X det 

fee7o 



a 6 70 



aGao 



/^(Aa, Aj) 
h{\j,\a) 



X det 

fce7Q 



a£ao 



■ (C.7) 



It is easy to see that this sum over partitions can be written as a single multiple integral 



Yi^Yl ^^^^^^ ^^^i'^'^^^j^ 



Ni f ,^42^^11 d{\,)d{^i,) 

r(Mu{A}) J-^ 



detjv 0^(14, MK^) detjv n^{{z}, {\}\{z}) ^ ^^ ^^ 

Here the integration contour r({/i} U {A}) surrounds the poles of $7^ in the points {A} and {;u} 
and no other pole of the integrand. 

References 

[1] H. Bethe, Zeitschrift fiir Physik, 71 (1931) 205. 

[2] R. Orbach, Phys. Rev., 112 (1958) 309. 

[3] L. R. Walker, Phys. Rev., 116 (1959) 1089. 

[4] M. Gaudin, La Fonction d'Onde de Bethe, Paris: Masson, 1983. 

[5] E. Lieb, T. Shultz and D. Mattis, Ann. Phys., 16 (1961) 407. 

[6] E. Lieb and D. Mattis (eds.). Mathematical Physics in One Dimension, New York: Aca- 
demic Press, (1966). 

[7] C. N, Yang and C. P. Yang, Phys. Rev., 150 (1966) 321, 327. 



19 



[8] L. D. Faddeev, E. K. Sklyanin and L. A. Takhtajan, Theor. Math. Phys. 40 (1980) 

[9] V. E. Korepin, N. M. Bogoliubov, A. G. Izergin, Quantum Inverse Scattering Method and 
Correlation Functions, Cambridge University Press, 1993. 

[10] B. M. McCoy, Phys. Rev., 173 (1968) 531. 

[11] A. Lenard, J. Math. Phys. 7 (1966) 1268. 

[12] M. Sato, T. Miwa, M. Jimbo, Publ. Res. Int. Math. Sci. 14 (1978) 223; 15 (1979) 201, 577, 
871; 16 (1980) 531. 

[13] V. E. Korepin and N. A. Slavnov, Comm. Math. Phys. 129 (1990) 103. 

[14] F. Colomo, A. G. Izergin and V. E. Korepin and V. Tognetti, Theor. Math. Phys. 94 (1993) 
11. 

[15] A. G. Izergin, V. E. Korepin, Comm. Math. Phys. 94 (1984) 67 

[16] A. G. Izergin and V. E. Korepin, Commun. Math. Phys. 99 (1985) 271. 

[17] V. E. Korepin, Comm. Math. Phys. 113 (1987) 177. 

[18] T. Kojima, V. E. Korepin and N. A. Slavnov, Comm. Math. Phys. 189 (1997) 709. 

[19] N. Kitanine, J. M. Maillet and V. Terras, Nucl. Phys. B, 554 [FS] (1999) 647, 



math-p h/9807020| 



[20] N. Kitanine, J. M. Maihet, and V. Terras, Nucl. Phys. B 567, 554 (2000), 
[^ath-ph/9907019 . 

[21] N. Kitanine, J. M. Maillet, N. A. Slavnov, V. Terras, in "Solvable lattice models 2004", 
RIMS pubhcations, 2006, ,hep-th/0505006, 



[22] J. M. Maillet and V. Terras, Nucl. Phys. B 575, 627 (2000), |hep-th/99H030| 



[23] N. Kitanine, J. M. Maillet, N. A. Slavnov, V. Terras, Nucl. Phys. B712 (2005) 600, 
|hep-th/0406190[ 

[24] V. E. Korepin, Comm. Math. Phys. 86 (1982) 391. 

[25] A. G. Izergin, Sov. Phys. Dokl., 32 (1987) 878 

[26] N. A. Slavnov, Theor. Math. Phis. 79 (1989) 502. 

[27] N. A. Slavnov, Zap. Nauchn. Sem. POMI 245 (1997) 270. 



20 



[28] N. Kitanine, J. M. Maillet, N. A. Slavnov, V. Terras, Nucl. Phys. B 641 [FS] (2002) 487; 
|hep-th/0201045 [ 

[29] N. Kitanine, J. M. Maillet, N. A. Slavnov, V. Terras, Nucl. Phys. B729 (2005) 558; 
[hep-th/0407108 , 



21 



